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ABSTRACT 


Subsurface radar, also known as ground-penetrating radar is increasingly being 
used for the detection and location of buried objects such as mines and structures that are 
found within the upper regions of the earth’s surfaee. The thesis gives a review of the 
work done to date in this area, laying emphasis on the possible antenna designs to match 
the range of intended applications. An overall design strategy is outlined, together with a 
more detailed treatment of the ground-penetrating radar subsystems and topics which are 
relevant to effeetive subsurface radar operation. These include the dieleetric properties of 
earth materials, the ehoiee of the frequeney of operation, as well as the design and con¬ 
struction of suitable antennas. 

Finally, a new antenna structure ealled the counter-wound spiral antenna, whieh is 
suitable for subsurface radar applications, is examined. The counter-wound spiral antenna 
has a broad bandwidth and a linear polarization with a controllable plane of polarization 
from a planar geometry. It has an electronieally steerable plane of polarization. This 
unique property offers a reduction in antenna polarization loss and allows the extraetion 
of maximum information from the target scattered eeho. 
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EXECUTIVE SUMMARY 


Subsurface radar or ground-penetrating radar (GPR) refer to a wide range of elee- 
tromagnetie teehniques designed primarily for the loeation of objeets or interfaces buried 
beneath the earth’s surfaee. These teehniques offer rapid, high resolution and non- 
invasive investigation of underground objeets and struetures by reeording mierowave ra¬ 
diation that passes through the ground, refleets from the buried objects or interfaces and 
is then returned to the surfaee. 

The operational effeetiveness of the subsurfaee radar depends on the (z) effieient 
eoupling of eleetromagnetie radiation into the ground; (zz) adequate penetration of the ra¬ 
diation through the ground having regard to target depth; {in) ability to obtain from bur¬ 
ied objeets, or other dielectrie discontinuities, a sufficiently large scattered signal for de- 
teetion at or above the ground surfaee; and (z'v) an adequate signal bandwidth in the de- 
teeted signal with regard to the desired resolution and noise levels. 

In GPR applieations, very often the antenna is placed directly over the ground 
with its axis normal to the surfaee without eontaeting the surface of the earth. This allows 
rapid surveying of the subsurface during operation. For this applieation, antennas with 
uniform properties such as input impedance, gain, and polarization over a wide band¬ 
width are desirable, espeeially when the transmitted and reeeived signal has a very large 
signal bandwidth. 

It has been observed experimentally that the attenuation of eleetromagnetie waves 
inereases with frequeney and that for a partieular frequeney, wet materials exhibit higher 
attenuation as eompared to dry ones. Variations in eleetrieal properties of ground materi¬ 
als eause refleetions. It is the refleetions from buried objeets that propagate baek to the 
surfaee that are reeorded by the subsurfaee radar. 

The requirements of the antenna designed for subsurface radar applieation ean be 
quite different from those that are used for atmospherie teleeommunieation purposes. The 
subsurfaee radar antenna has to be matehed to the eharaeteristies of the propagation me¬ 
dium as well as the geometrieal shape of the target. 
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The impulse response, fractional bandwidth and polarization state of the antenna 
are some of the critical parameters being considered. The cross-coupling levels between 
the closely spaced transmitter and receiver, the interaction between the reactive field of 
the antenna and the propagation medium as well as the geometry (as in planar or non- 
planar) of the antenna are also important considerations. 

The spiral antenna, which belongs to the class of frequency-independent antennas, 
has radiation pattern, impedance and polarization that remain virtually unchanged over a 
large bandwidth. An ideal spiral antenna that radiates a circularly polarized wave towards 
the air/ground interface will not receive the wave reflected from the interface, because it 
is predominantly cross-polarized. When used for detection of objects such as a rod that 
scatters a linearly polarized wave, the spiral antenna only receives the co-polarized half 
of the backscattered power, therefore reducing the probability of target detection. 

For counter-wound spiral antenna (CWSA), the polarization of the wave transmit¬ 
ted by the antenna can be made either vertically or horizontally polarized by varying the 
delays to the signals that are being sent to the antenna feeds of the spirals. This design 
allows buried objects to be detected by both vertically and horizontally polarized waves, 
which results in reduced polarization loss. This enables better detection and identification 
of targets which lie in any arbitrary orientation with respect to the main antenna axis. 

Several important aspects of the performance of the CWSA have been verified 
through simulations. It has been shown that the rotated CWSA (with the plane of the 
front and rear spiral antenna planes rotated 90 degrees to each other along the main an¬ 
tenna axis), gives a better performance in terms of reflections as compared to the un¬ 
rotated case. The CWSA far field polarization has also been verified to be linearly polar¬ 
ized when an appropriate delay is set to the excitation input of the front planar antenna. 
The plane of polarization of the linearly polarized wave can be changed by varying the 
signal excitation delay. The research work has resulted in the successful design of a 
CWSA that has a broad bandwidth and linear polarization with an electronically control¬ 
lable plane of polarization from a planar geometry. 
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I. INTRODUCTION 


The terms subsurface radar or ground-penetrating radar (GPR) refers to a wide 
range of electromagnetic techniques designed primarily for the location of objects or in¬ 
terfaces buried beneath the earth’s surface. These techniques offer rapid, high resolution 
and non-invasive investigation of underground objects and structures by recording mi¬ 
crowave radiation that passes through the ground and is returned to the surface. Virtually 
all systems developed in recent years to detect landmines and buried unexploded ordi¬ 
nance (UXO) have included ground-penetrating radar (GPR) as one of the principle sub¬ 
systems. GPR is currently one of the leading sensor technologies the U.S. Army is inves¬ 
tigating for the detection of sub-surface landmines. Planning System Inc, has developed 
working stepped frequency GPR systems that operate over the frequency band of 800 
MHz to 4 GHz, using different variations of log spiral and Archimedean spiral antennas 
[ 1 ]. 

In a typical GPR, a transmitter sends a microwave signal into the subsurface 
which will be reflected by the buried objects or interfaces beneath the earth’s surface. 
The microwaves propagate at velocities that are dependent upon the dielectric constant of 
the subsurface medium. Changes in the dielectric constant that are due to changes in the 
subsurface materials cause the radar wave to be reflected. The time taken for the energy 
to return to the surface is related to the depth at which the energy was reflected. Thus in¬ 
terpretation of this reflected energy yields information on structural variation of the sub¬ 
surface. 

The design of GPR is largely applications-oriented and the overall design phi¬ 
losophy is usually dependent on the target type and the obscuring medium. For example, 
the choice of GPR operating frequency is correlated to the expected depth of penetration 
of the microwave signals and the required ground vertical resolution. Higher frequency 
sources will offer greater vertical resolution of structure, but will not penetrate as deep as 
lower frequency sources. 

The range of applications of subsurface radar methods is growing, along with the 
sophistication of signal recovery techniques, hardware designs and operating practices. 
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with the last including airborne and satellite surveying. GPR is the only geophysical tech¬ 
nique that can offer the horizontal and vertical resolution necessary for many ground 
probing applications. This includes hydrological and geotechnical applications such as 
determining soil and bedrock characteristics in the shallow interface, as well as high- 
resolution imaging for monitoring structure integrity of buildings, mine walls, roadways 
and bridges. 

The first description of the use of electromagnetic signals to determine the loca¬ 
tion of buried objects was found in a German patent by Leimbach and Lowy [2] in 1910. 
Their technique consisted of burying dipole antennas in an array of vertical boreholes and 
comparing the magnitude of signals received when successive pairs were used to transmit 
and receive. The main feature of this technique is the use of continuous wave (CW) to¬ 
gether with diffraction effects of buried objects or underground features. Scattering from 
objects are due to variation in the conductivity of the subsurface medium. 

The first use of pulsed techniques to determine the structure of buried features ap¬ 
pears to originate from the work of Hulsenbeck in 1926 [3]. Pulsed techniques were then 
extensively developed over the years as a means of probing to considerable depths in 
subsurface media including ice, fresh water, salt deposits, desert sand, and rock forma¬ 
tions. A more extended account of the history of subsurface radar and its growth up to 
about 1976 is given by Daniels [4] and Nilsson [5]. 

Alternative techniques that have been used for detecting buried objects or under¬ 
ground structures include the seismic method and the use of sonar in marine environ¬ 
ments. Seismic methods are used extensively for the exploration of deep geographical 
features. Other methods include resistivity and magnetometry surveying, as well as low- 
frequency induction for detecting localized shallow metal objects or deeper metallic pipes 
and cables. The preferred method to be used depends on the target type, the nature of the 
immersing material as well as the operational requirements. For detection of non-metallic 
targets, a seismic method is a good alternative to subsurface radar. 

One of the operational advantages of using subsurface radar techniques is its abil¬ 
ity to perform rapid surveying of the subsurface. This is due to the absence of contact be- 
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tween the antenna and the surfaee of the earth. The GPR antenna ean also be designed to 
achieve adequate bandwidth to suit its intended application. 


In a GPR applications, very often the antenna is placed directly over the ground 
with its axis normal to the surface without contacting the surface of the earth. This allows 
rapid surveying of the subsurface during operation. For this application, antennas with 
uniform properties such as input impedance, gain, and polarization over a wide band¬ 
width are desirable, especially when the transmitted and received signal has a very large 
bandwidth. The two-arm equiangular spiral antenna discussed here has these frequency- 
independent characteristics in free space, and these characteristics are preserved to a great 
degree when placed over the ground. Olhoeft [6] discussed several applications and frus¬ 
trations in using ground-penetrating radar system. 

This thesis presents a review of subsurface radar, highlighting those features 
which differ from conventional atmospheric radar techniques and emphasizing the inter¬ 
relationship between those topics in electromagnetism, soil science, geophysics and sig¬ 
nal processing, which form a critical part of the design of a GPR system. A detailed dis¬ 
cussion of the design of a two-arm frequency-independent equiangular spiral antenna is 
presented, showing the important antenna characteristics. A simple design of a 4-arm spi¬ 
ral antenna and its characteristics is then discussed. 

The thesis is organized as follows. Chapter II gives a review of the work done to 
date in the area of subsurface radar, the overall design strategy together with a more de¬ 
tailed treatment of the ground-penetrating radar subsystems and topics which are relevant 
to effective subsurface radar operation. These include the dielectric properties of earth 
materials, signal modulation schemes and design and construction of suitable antennas. 

Chapter III gives a detailed discussion of the construction of a dual-arm fre¬ 
quency-independent planar spiral antenna, comparing its performance with the theory. 
Chapter IV presents the design of counter-wound spiral antenna, illustrating its perform¬ 
ance with simulation results from Microwave Studio. Chapter V presents a summary and 
proposes prospective developmental work in this arena. 
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II. OVERVIEW OF SUBSURFACE RADAR TECHNIQUES 


This chapter examines the principal features affeeting the operation and design of 
a ground-penetrating radar. The deseriptions here are mainly qualitative. More detail de¬ 
sign information can be found in the teehnical papers refereneed. 

The monostatic radar technique is usually employed to deteet back-scattered ra¬ 
diation or refleetion from a buried target. Although forward scattering of the target ean 
also be used, it is seldom employed due to the need for at least one antenna to be buried 
and some kind of image transformation to be applied to the measured data. All of the 
work diseussed here is based on target backscattering measurement. 

Figure 1 shows the block diagram of a general subsurface radar system. The de¬ 
sign is typieal of a base-band, sub-surfaee pulsed radar system. The subsurfaee radar sys¬ 
tem reeords mierowave radiation that passes through the ground and is returned to the 
surfaee. A transmitter sends a mierowave signal into the subsurfaee. As the radar wave 
propagates at a veloeity that is dependent upon the dielectric constant of the subsurface 
medium, changes in the dieleetric eonstant that are due to changes in the subsurface ma¬ 
terials will eause the radar wave to reflect. The time it takes for the electromagnetie en¬ 
ergy to return to the surfaee is related to the depth at which the energy was reflected. In¬ 
terpretation of this reflected energy yields information on structural variation of the near 
subsurface or the depth of the buried objects. 

Separate broadband antennas are usually used for transmission and reeeption. The 
design options available are usually dietated by the electromagnetie properties of the 
ground, while the target type and shape will affeet the choice of antenna type and eon- 
figuration as well as the kind of proeessing teehniques to be employed. The effects of the 
ground on the antenna can be largely deeoupled from the rest of the problem and solved 
independently. 

The operational effectiveness of the subsurfaee radar depends on (z) the efficient 
coupling of electromagnetic radiation into the ground; (zz) adequate penetration of the ra¬ 
diation through the ground having regard to target depth; (zzz) the ability to obtain from 
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buried objects, or other dielectric discontinuities, a sufficiently large scattered signal for 
detection at or above the ground surface; and (iv) an adequate signal bandwidth in the 
detected signal with regard to the desired resolution and noise levels. The propagation 
loss, clutter characteristics and target characteristics described here are distinctly different 
from those of the conventional free-space radar system [4], and the tradeoffs involved are 
significantly different. 



Figure 1. Block diagram showing the operation of a typical base-band subsurface 

radar system (From Ref. [4].). 

The quasi-monostatic geometry of the subsurface radar is shown in Figure 2. The 
transmitting and receiving antenna are located at almost the same elevation above the 
ground at which the object/scatterer is buried. The antenna transmits an incident field to¬ 
wards the ground. This field is partially transmitted into the ground and partially reflected 
from the air/ground interface. The wave that is transmitted into the ground is then scat¬ 
tered by the buried object. In the figure, represents the radar cross section (RCS) of 
the target, R is the slant range, d is the depth of the buried object, Tx and Rx are the 
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transmitting and receiving antennas respectively, and 0 is the incident angle of the trans¬ 
mitted wave measured from the z-axis. 


Tx Rx 



Figure 2. Quasi-monostatic geometry of a typical subsurface radar system 

(From Ref [7].). 


A, DIELECTRIC PROPERTIES OF EARTH MATERIALS 

Ground material can be described by its electrical properties such as electric per¬ 
mittivity (s), electric conductivity (cr) and magnetic permeability (ju). Since most of 
the earth’s material is non-magnetic, the permeability of the medium can be approxi¬ 
mated as the permeability of free space {ju = jUq). The conductivity of the medium deter¬ 
mines the amount of energy lost from currents induced in the medium by the propagation 
wave. 

The ground materials are generally complex, and losses can arise from both con¬ 
duction and displacement currents. In the first case, ohmic loss is conventionally ex¬ 
pressed in terms of the material conductivity, whereas for a dielectric, a complex permit¬ 
tivity is used. When both losses are present (for example, in wet sand) the losses are not 
kept separate, but lumped together in either an equivalent conductivity or equivalent 
permittivity. 
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Therefore several different representations are possible. Complex permittivity (^) 


and complex conductivity (cr) can be defined as follows 


complex permittivity 

s = s'-js" 

(la) 

complex conductivity 

f • ft 

(j = G -J<J 

(lb) 

effective permittivity 

II 

r 

(2a) 

effective conductivity 

<j'^=(j' + cos” 

(2b) 

imaginary effective permittivity 

05 

II 

+' 

(2c) 

imaginary effective conductivity 

cr" = cr" - 1 - cos'. 

(2d) 


When measurements are made on a conducting dielectric, the parameters meas¬ 
ured are the apparent permittivity (?) and the apparent conductivity (a), defined as 




(3a) 


? = o-; 




(3b) 


The effective electric loss tangent (tan?^) is given by 


tanS^ = 


cr' + cos” 
cos' -a" 




(4) 


The measurements of interest for making an assessment of a material type for 
subsurface radar applications are its attenuation («) and the velocity of wave propaga¬ 
tion (o) as a function of frequency. The relationship between these parameters and the 
apparent permittivity (?) is given by the equation 

r = a + j/3 = . (5) 

The dielectric constant defines the index of refraction of the medium and is the 
material constant which controls the phase velocity of the electromagnetic waves in the 
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material. When the permeability is taken to be the permeability of free spaee (//^), the 
wave velocity is given by 


o = c 


2s, 


■('\/l+~tan^~^ + 1) 


( 6 ) 


where c is the speed of light in air. From (2), changes in the subsurface material will af¬ 
fect the index of refraction, and reflected energy will be produced related to the contrast 
in the dielectric constant across the boundary between the two materials. 

The attenuation constant is given by 





(V) 


or 


a[dB/m] = 129/^/^ ^(1 + tan^^J-l 


( 8 ) 


where / is in GHz ((y = 2;z-/). 

Table 1 shows the dielectric constant of some common earth materials. The di¬ 
electric constant of the medium depends principally on the water content. The reason is 
that at lower frequencies, water has a relative permittivity of about 80, while for the solid 
constituents of most soils, s^ is in the range of 2 to 6. The measured s^ of soil was found 
to be in the region of 4 to 40 [1]. 
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Material 

Dielectric Constant 

Propagation Velocity 

(m/ns) 

Air 

1 

0.3 

Water 

81 

0.033 

Granite 

9 

0.10 

Limestone 

6 

0.12 

Sandstone 

4 

0.15 

Rocks 

4~12 

0.15-0.087 

Dry sand 

4~6 

0.15-0.12 

Wet sand 

30 

0.055 

Dry clay 

8 

0.11 

Wet clay 

33 

0.052 

Dry soils 

3~8 

0.17-0.11 

Wet soils 

4~40 

0.15-0.047 

Asphalt 

3~6 

0.17-0.12 

Concrete 

9~12 

0.10-0.087 


Table 1. Typical dielectric constant of common soil medium and velocity of wave 

propagation (From Ref [I].)- 

It has been observed experimentally that, for most materials that make up the 
earth near its surface, the attenuation of the electromagnetic radiation increases with fre¬ 
quency. In general, wet materials exhibit higher loss than dry ones at a given frequency. 
Figure 3 below shows the one-way path loss for different material surfaces over a fre¬ 
quency range from 1 MHz to 1 GHz. It is shown here that the attenuation of wet rock is 
higher than dry rock, and salt water has the highest level of attenuation. 
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Figure 3. One-way path loss for different surface materials over 1 GHz frequency 

range (After Ref. [7].). 


There is sufficient information to enable a realistic assessment to be made of the 
likely range of dielectric properties to be encountered when a subsurface radar is being 
designed for work in a particular environment. Table 2 gives an account of published 
work relating either to the experimental or the theoretical studies of the dielectric nature 
of a range of material. 

Variations in electrical properties cause reflections at high frequency. The theory 
of how the reflected waves and the transmitted (refracted) waves propagate away from 
the boundary of different media is given by Snell’s Law. It is the reflections from buried 
objects that propagate back to the surface that are recorded by the subsurface radar. 
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Material 


Research Area /References 


Water, lee, snow 

Properties of water in solid and liquid state, King & Smith [8], 


Frequency dependence of water dielectric properties, Potter [9], 


The complex-dielectric constant of sea-ice at frequencies, Vant [10], 


Snow properties effects on millimeter-wave backscatter, Larry [11], 


Loss-Factor behavior of freshwater ice, Georgi [12], 



Rook, minerals 

Dielectric properties of rocks, King & Smith [8]. 


Data on dielectric properties of rocks, Parkhomenko [13], Keller [14], 

Ulaby [15], 


Radio Propagation at 900 MFlz in underground coal mines, Zhang [16], 


Dielectric properties of salt deposits, Unterberger [17] 



Soils 

Low frequency (<100 Flz) soil dielectric characterization, Dmitriev[18], 

Fano [19], 


Moisture Effects on the dielectric properties of soils, John [20]. 


Dielectric Properties of soils in the 0.3-1.3 GFlz , Neil [21]. 


Complex dielectric constant of wet and frozen soil, Tikkhonov [22]. 




Table 2. Published work relating to the experimental and theoretieal studies of the 

dieleetrie properties of a range of materials. 
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B. DEPTH AND FREQUENCY OF OPERATION 

Subsurface radar systems usually operate in the megahertz range. The waves that 
propagate into the ground will have wavelengths on the order of 1 meter or less. The 
horizontal and vertieal resolutions of the radar images are dependent upon the wavelength 
of operation, sueh that the smaller the wavelength, the better the resolution. The higher 
frequeney souree will, therefore, yield better resolution, but the higher frequeney signals 
will not penetrate as deep as the lower frequeneies. For a given signal deteetion threshold 
the maximum depth of investigation deereases rapidly with inereasing frequeney. Thus a 
eareful ehoiee of frequeney of operation must be made based on the expeeted target depth 
and the system design goals. 

Table 3 shows the typieal required penetration depths of a range of materials, and 
their appropriate upper operating frequencies. A eonsideration whieh applies when 
ehoosing an operating frequeney for deteeting loealized objeets, rather than interfaees 
between thiek layers, is the faet that the baekseattered amplitude is also frequeney de¬ 
pendent and, apart from any resonanee effeets, deereases with deereasing frequeney. 

It ean be eoneluded from Table 3 that subsurfaee radar operating in the 100 MHz 
to 1 GHz frequeney range has a penetration depth of about 2 to 3 meters into the different 
types of soils (elay, loam and sand). This will be the typieal operating frequeney range of 
most subsurfaee radars that are used for mine deteetion. 

In longer range (deeper) investigation of geologieal strueture or deep iee, frequen¬ 
eies as low as a few tens of MHz are used. This is as shown in the ease of 20 meters of 
rook or ooal, with the radar operating at a muoh lower frequeney in the 50 MHz to 500 
MHz range. 

The deepest penetration therefore ooours in dry, non-clay soils, and in dry rocks 
with no clay cementation. Snow and iee oover will have little effeets on GPR data oolleo- 
tion. When the soils or rooks are saturated, the oonduotive nature of the filling liquid be- 
oomes important. A eonsideration when ehoosing an operating frequeney for the deteo- 
tion of loealized objeets rather than interfaees is that the baoksoatter from the targeted 
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object is also frequency dependent. The amplitude of the backscatter will also decrease 
with increased frequency of operation. 


Material 

Penetration depth 

Max. Operating Freq 

Cold pure fresh water ice 

10 km 

10 MHz 

Temperate pure ice 

1 km 

2 MHz 

Saline ice 

10m 

100 MHz 

Fresh Water 

100 m 

100 MHz 

Sand (desert) 

5 m 

1 GHz 

Sandy soil 

3 m 

1 GHz 

Loam soil 

3m 

500 MHz 

Clay soil 

2 m 

100 MHz 

Salt (dry) 

1 km 

250 MHz 

Coal 

20 m 

500 MHz 

Rocks 

20 m 

50 MHz 

Walls 

0.3 m 

10 GHz 


Table 3. Depths at which radar probing gives useful information, taking into account the 
attenuation encountered and the nature of the reflectors (From Ref [4].). 

C. DEPTH RESOLUTION 

The depth resolution required for the subsurface radar is dependent on its intended 
application. In situations where only the depth of the interface is of interest, the radar 
may only need to operate at a single frequency. The depth of the interface can be deter¬ 
mined by measuring the time delay of the returned echoes. This technique is often known 
as radio echo sounding. 


14 






For detection of specific buried objects such as mines or buried pipes, a larger 
bandwidth is required in order to be able to distinguish between the various scatterers. 
The earth, which acts as a low pass filter, would modify the spectrum of the transmitted 
signal according to its material properties. This effect of a frequency-dependent attenua¬ 
tion of the electromagnetic waves by the ground limits the information bandwidth of the 
received signal at a given noise or clutter level. 

Figure 4 shows the spectral spread of a stepped frequency waveform (12 pulses, 
frequency step size 50 MHz per pulse, start frequency 500 MHz, PRF 10 MHz, pulse 
width 0.01 microsecond) being transmitted and received from a target buried under¬ 
ground. 



Figure 4. Frequency spectral spread of a stepped frequency waveform being 
transmitted into and received from the ground (After Ref. [7].). 
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From Figure 4, it can be observed that the ground loss has reduced the effective 
bandwidthby about 25%, where is given by 

^00 ^ 

S/=-j(2Vr|S(/)|#, (9a) 

—00 

where E is the signal energy and the range resolution of the radar is approximated by 

(9b) 

where is the effective pulse width. The ground loss has also caused an increase in the 
range cell by about 25%. 


D, HORIZONTAL RESOLUTION 

The horizontal resolution of the subsurface radar is important when localized tar¬ 
gets of the same depth are being sought and distinguished. The technique employed to 
obtain horizontal resolution depends on the wave attenuation of the medium. The loca¬ 
tions of the transmitter and receiver, as to whether they are collocated or physically sepa¬ 
rated, will have to be known. 

The horizontal resolution of the subsurface radar can be found in several ways. 
One is based on the power distribution of the received signals from the scatterers as the 
radar moves over the buried target. Figure 5 shows a collocated transmitter and receiver 
pair being scanned along the ground surface above a buried object at depth d, with the 
one way voltage attenuation constant of the medium as a . The monostatic case is shown 
with the incidence angle of the transmitted wave measured from the z-axis denoted as 6. 
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Tx/Rx 



Figure 5. Geometry for estimation of the horizontal resolution of subsurface radar 

(From Ref [7].). 


The slant range R can be calculated if the depth d of the buried object and the 
distance x from the z-axis, directly above the buried object is known. Range R is given 
by 

R = . (9) 

The field (voltage) attenuation constant of the medium a is 

a = -ln(10'^'''’), (10) 


where L is the “ loss in dB per meter” of the medium, and is related to the attenuation 
constant by the equation 


L = 20 logf = 20 logf■). 


E at surface 


( 11 ) 


Assuming an inverse fourth power of decay of received signal power with target depth 
(i.e., the monostatic case), the signal received from the target varies as 

-4 off 

(12) 
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Figure 6 shows the return from a target at depth d = 2 m, with the radar moving 
on the ground along the x-axis. It shows the normalized power reeeived from the target 
for three ground attenuation values. Figure 7 shows the returns from two equal RCS 
targets positioned at x = 0 and x = 2 m, both at depth d = 2 meters. The normalized re¬ 
ceived powers from different position of x for two different ground losses of 3 dB/m and 
10 dB/m were simulated. It was observed that higher ground loss would actually im¬ 
prove the horizontal resolution. Therefore, increasing the loss per meter of the medium 
has the same effect as narrowing the antenna beamwidth, resulting in better horizontal 
resolution. 



Figure 6. Power received from target at depth {d=2vs\) for varying values of 

distance x away from the target in three types of ground losses (After Ref. [7].). 
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Solid line : ground loss of 3 dB/m 
Dashed line : ground loss of 10 dB/m 



Figure 7. Power received from 2 target at depth (t/ = 2 m) and separated by x = 2 m, 
for ground loss of 3 dB/m and 10 dB/m (After Ref [7].)- 
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E, FREQUENCY-INDEPENDENT ANTENNAS 

The requirements of an antenna designed for subsurface radar application can be 
quite different from those that are used for atmospheric telecommunication purposes. The 
subsurface radar antenna has to be matched to the characteristics of the propagation me¬ 
dium as well as the geometrical shape of the target. 

The impulse response, fractional bandwidth, and polarization state of the antennas 
are some of the critical parameters being considered. The cross-coupling levels between 
the closely spaced transmitter and receiver, the interaction between the reactive field of 
the antenna with the propagation medium as well as the geometry (as in planar or non- 
planar) of the antenna are also important considerations. 

As attenuation of radio waves in geophysical media (such as soil or ice) increases 
with frequency [4, 23], better ground penetration is usually achieved by operating the ra¬ 
dar at lower frequency. This would however results in bigger size antenna since antenna 
size are generally proportional to the frequency of operation. 

The range resolution of the radar depends on its operating bandwidth. The wider 
the frequency range of operation, the better would be its range resolution. Thus, to 
achieve finer range resolution, a broadband antenna is required. Typical ground- 
penetrating radar with a 500-MHz frequency range of operation would correspond to a 
radar range resolution of about 15 cm in soil [23]. The lowest frequency of operation is 
usually a compromise between the media attenuation and the antenna size [24]. 

Element antennas like monopoles, cylindrical dipoles, bi-conical dipoles and re- 
sistively loaded bowtie antennas are characterized by limited bandwidths, low directivity 
and linear polarization. The antennas are, however, usually small in size, and can be ar¬ 
ranged cross-polarized to provide good isolation and effective target discrimination 
against planar targets. The current distribution of the element antennas can be conven¬ 
iently approximated by current elements, making the analysis of the radiation characteris¬ 
tics convenient. 

A circularly polarized antenna is usually preferred to a linearly polarized antenna 
because the strength of the reflected wave from the object does not depend on the orienta- 
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tion of the object relative to the antenna. For a linearly polarized antenna, if the transmit¬ 
ting antenna is cross-polarized to the receiving antenna in order to reduce mutual cou¬ 
pling, the receiving antenna will probably not be able to detect the reflected wave from 
thin object. 

Frequency-independent antennas like the Achimedian (or sometimes known as 
Archimedes Spiral), equi-angular spiral and conical spiral are all circularly polarized 
broadband antennas. Dyson [25] discussed the basic differences between the equiangular 
and Archimedes spiral antenna and their characteristics, such as radiation efficiency and 
ultimate bandwidth capabilities. 

Among the circularly polarized antennas mentioned above, only the conical spiral 
antenna radiates unidirectionally. However, on the conical spiral antenna, different fre¬ 
quency waves are radiated from different regions of the cone. This gives rise to problems 
in range measurement, because the distance to the target changes with frequency. Hertel 
[26] investigated the two-arm conical spiral antenna when it is placed directly over the 
ground with its axis normal to the surface of the ground. The conical spiral antenna has 
been shown to have uniform impedance, gain, and circular polarization over a broad fre¬ 
quency range. 

Planar antennas like the equi-angular spiral, Achimedian, and log-periodic spiral 
antennas, shown in Figures 8 through 10 do not have the range measurement problem 
mentioned earlier. However, the planar antenna is bi-directional. Takashi [24] demon¬ 
strated the use of planar spiral antenna with a cylindrical waveguide back to achieve a 
wideband (500 MHz bandwidth), circular polarization and unidirectional properties. The 
minimal physical volume of planar antennas is certainly a distinct advantage. 
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Figure 8. Four-arm logarithmic spiral antenna with self-complementary pattern 

(From Ref. [4].). 



Figure 9. Two-arm Archimedean spiral antenna (From Ref. [4].). 



Figure 10. Four-arm planar log-periodic antenna (From Ref [4].). 
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Subsurface radar usually employs separate antennas for transmission and reeep- 
tion. This is due to the absenee of fast switehes to protect the receiver during transmis¬ 
sion. The cross-eoupling level between the transmitting and reeeiving antennas (isolation) 
is an important eonsideration. 

Some of the general design principles of frequeney-independent antennas dis- 
eussed by Rumsey [27] are as follows: 

a) The antenna is excited from the end at which the high frequencies eleetro- 
magnetic waves are radiated. 

b) The inaetive part of the antenna should produce negligible radiation. 

c) The geometry of the antenna should be defined entirely by angles. 

Many of the antennas are fed by a balun, to match the unbalaneed transmission line to a 
balanced radiating strueture. 

Other types of antennas that have been successfully used for subsurface radar ap¬ 
plications are traveling-wave antennas and aperture antennas. The classes of antenna that 
can be used successfully for subsurfaee radar applications generally require a wide band¬ 
width and low operating frequeney. The choice between a planar and non-planar geome¬ 
try is more eritieal when the antenna is required to operate elose to the ground surface. 
The equiangular spiral antenna is useful in eases where an extended impulse response ean 
be aeeepted but, under conditions when the impulse response must be more limited, the 
eontinuously-loaded dipole or traveling-wave antenna is preferred. 


F. DATA COLLECTION 

The most common type of survey technique used in subsurface radar applieation 
is common offset profiling. Other techniques inelude the common midpoint survey. In the 
eommon offset profiling teehnique, the transmitting and reeeiving antennas are at a fixed 
distance and the radar is progressively moved along the line of survey to reeord baekseat- 
tered returned signals from the subsurfaee. With bistatie antennas, the antenna separation 
(souree-reeeiver offset) has to be held eonstant for common offset profiling, and this off¬ 
set ean be optimized for best results. For this type of system, the horizontal sample inter- 
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val is usually chosen to be large whieh results in redueed lateral resolution. With a 
monostatie or quasi-monostatie system, both the transmitting and reeeiving antennas are 
housed within the same unit. This allows for better horizontal resolution profiling of the 
subsurface. 

In eommon midway profiling, the two bistatic antennas are progressively moved 
away from eaeh other, keeping the eenter between the two antennas fixed a one location. 
As the antennas are moved away from the origin, the subsurfaee data is eollected at eaeh 
new, more distant position. Sinee the eleetromagnetie wave travel time is direetly pro¬ 
portional to the inereasing offset between the two antennas, this type of data profiling can 
help to determine the velocity of the wave in the subsurface. 

During data aequisition, multiple data seans are usually taken at eaeh survey loea- 
tion. The data reeorded are then added up coherently to remove ineoherent noise present 
in the data samples. The monostatic mode discussed is usually used for rapid surveying in 
a continuous data profiling mode, whereas the bistatic mode is used for enhaneed resolu¬ 
tion profiling over discrete locations. 


G. DATA PROCESSING 

Some basic data processing funetions are usually performed on the raw data eol- 
leeted. The time offset due to instrument reeording is first removed prior to any data in¬ 
terpretation of the radar images. The ringing due to the close proximity of the souree and 
reeeiving antennas in monostatic system are then removed. The effects due to non¬ 
constant motion along the profile and amplitude variation would have to be compensated. 
The return signal due to reflections from near objects will naturally have a much higher 
amplitude as eompared to the those objects buried at a much deeper depth. This loss of 
signal amplitude is due to both the intrinsie attenuation of the soil and the geometrie 
spreading of the signals. Various time-variable gain funetions may be applied to eompen- 
sate for sueh gain variations in the reeeived signals. Lateral averaging of data by sum¬ 
ming data traces reeeived from the subsurfaee can be used to improve the signal-to-noise 
ratio of the reeeived signals. Frequency filtering is also used to remove unwanted spuri¬ 
ous signals present in the data samples to produce a more interpretable GPR image. 
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H. SUMMARY 


This chapter examined the principal features affeeting the operation and design of 
the ground-penetrating radar and the monostatic radar technique that is usually employed 
to deteet baek-scattered radiation or refleetion from a buried target. 

The operational effectiveness of the subsurface radar depends on several factors 
such as the efficient coupling and penetration of eleetromagnetic radiation into the 
ground and the ability to obtain a suffieiently large and wideband scattered signal. 

It has been observed that the attenuation of the electromagnetic radiation in¬ 
creases with frequency for most materials that make up the earth near its surface, and in 
general, wet materials exhibit higher loss than dry ones at a given frequency. 

The elasses of antenna that ean be used suecessfully for subsurface radar appliea- 
tions generally require a wide bandwidth and low operating frequency. The spiral antenna 
is a prime eandidate, and its performance will be examined in the next chapter. 
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III. SPIRAL ANTENNA DESIGN 


The design, characteristics, and performance of the planar spiral antenna is exam¬ 
ined in this chapter. The planar spiral antenna has a broad bandwidth and generates circu¬ 
larly polarized waves. These are useful properties for subsurface applications. 

A. PLANAR EQUIANGULAR SPIRAL ANTENNA IN FREE SPACE 

Frequency-independent antennas refer to antennas that have no theoretical limita¬ 
tion on the bandwidth of operation. In practice, these are antennas whose radiation pat¬ 
tern, impedance and polarization remain virtually unchanged over a large bandwidth [28]. 
Between the band limits, the performance varies in a manner that is periodic with the 
logarithm of frequency. These antennas are often also called as logarithmically periodic 
or log-periodic antennas. 

The design of a frequency-independent antenna is based upon the fundamental 
principle that, if the shape of the antenna were such that it could be specified entirely by 
angles, its performance would be independent of frequency [29]. 

This implies that an infinite electrical dimension of such antenna, in order to op¬ 
erate over an infinite frequency range. This means that the current distribution must de¬ 
crease with distance from the input much more rapidly so that the infinite antenna can be 
truncated with practically no effect on the radiation pattern [30]. To achieve frequency 
independence with a practical size antenna (one of finite size), the antenna must consist 
of active regions. An active region is that part of the antenna that contributes the most to 
the radiation at one particular frequency. Since this region moves along the antenna as the 
frequency is changed, its dimensions scale with wavelength and the antenna is able to 
fulfill the self-scaling and self-complementary conditions over a broad bandwidth. 

For the required frequency range of operation, the spiral antenna can be truncated 
while still retaining the properties of the infinite structures over a very wide bandwidth of 
operation. The arm length of the antenna need only be of the order of one wavelength at 
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the lowest frequency of operation for frequency independence. The upper frequency of 
operation is constraint by the size of the feed region [29]. 

Dielectric loading could be also be used to reduce the resonant frequency and the 
physical size of the antenna. This is done by placing a piece of dielectric substrate below 
the antenna [31]. In general, one is able to achieve better antenna efficiency and larger 
bandwidth by using a thick substrate with a low dielectric constant, at the expense of lar¬ 
ger element size. Similar performance can also be obtained using a thin structure that has 
a higher dielectric constant. 

Rumsey [27] has reviewed a wide range of frequency-independent antennas. The 
planar and conical equiangular spiral antennas have useful properties for subsurface ap¬ 
plications. Both produce circularly polarized radiation while four-element equiangular 
spirals can be designed to have two orthogonal pairs of elements, so that the use of one 
pair to transmit and one pair to receive allows the depolarization characteristics of targets 
such as pipes and tunnels to be exploited. 

The planar equiangular spiral antenna is an example of self-complementary struc¬ 
ture which enables its impedance to be calculated through Babinet’s principle [32]. For 
any self-complementary structure the antenna input impedance is frequency-independent 
in the designed bandwidth of operation. 

The planar equiangular spiral radiates in both directions perpendicular to the 
plane of the planar anteima. Absorbing material can be placed at the rear of the planar 
anteima if a unidirectional radiation pattern is required. Generally, it is not possible to 
obtain a wide bandwidth using a conducting back plane, because destructive interference 
between the spiral and its image will occur at some frequencies in the band. 

Thaysen [33] discusses the design of a cavity-backed coplanar-strip-fed logarith¬ 
mic single-plane spiral antenna that covers a 9; 1 bandwidth with a return loss better than 
10 dB operating from 0.4 GHz to 3.8 GHz. 

The electromagnetic simulation program Microwave Studio (developed by GST) 
was used to analyze the performance of the spiral antenna in terms of its input reflection 
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coefficients, antenna radiation pattern, and input impedance. The simulated results of the 
antenna were compared to theoretical results. 


B. PLANAR EQUIANGULAR SPIRAL ANTENNA GEOMETRY DESIGN 

The design of the equiangular spiral antenna is based upon a simple fundamental 
principle that if all dimensions of the antenna active regions are changed in linear propor¬ 
tion to a change in wavelength, the performance of the antenna is unchanged except for a 
ehange in scale in all measurements of length. The surface of the logarithmic planar spi¬ 
ral antenna ean be completely described by angles using the following equations 



(13a) 



(13b) 


where and Pj are the outer and inner boundaries of the spiral, respectively. The quan¬ 
tity Pq is the initial inner radius of the spiral, a is the growth rate, (j) is the angular posi¬ 
tion varying from 0°to 540° (for 1-1/2 turn spiral arm design), and is the initial angu¬ 
lar displacement between the two spiral arms. These dimensions are as shown in Figure 
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Figure 11. Geometry of the planar equiangular spiral antenna. 
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Since the spiral antenna was designed using the equation 




acf) 


(14) 


if the wavelength, X , is chosen as the unit of the antenna length and if p' equals the ra¬ 
dial coordinate measured in wavelengths, then 






(15) 


where 


^0 =-ln(;i). 
a 


(16) 


The above result shows that the effeet of ehanging wavelength is equivalent to changing 
the angle . The active region of the antenna in Figure 12 moves along the spiral to¬ 
wards the ends with decreasing frequency. Consequently, the highest frequency ( fmax) 
signal is radiated from the small end of the antenna and the lowest frequency {fmin ) from 
the large end. In the simplest approximation, the active region at a given frequency oc¬ 
curs when the spiral is roughly one wavelength in circumference. 

The bandwidth (BW) of broadband antennas is often defined as the ratio of the 
frequeneies that mark the bounds of acceptable performance, so the bandwidth of the 
equiangular spiral is given by 


BW = 


f 

J m 


f. 

J mm 


(17) 


Experimental and simulation results, however, show that the actual bandwidth is gener¬ 
ally lower than that obtained from Equation (17). The radiation pattern of the infinite 
structure of the spiral antenna is thus independent of frequency except for a rotation in 
the radiation pattern. This has also been found to be true for the finite length structure 
[29]. 

Thaysen [29] has determined a set of design parameters which will give the most 
frequency-independent radiation pattern and a fairly constant antenna input impedance. 
The parameters used for the geometry of a single planar spiral antenna were [29] 
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Po = 2.1mm 


(18a) 


a = 0.5 rad’' 


(18b) 


6»o =1.3 rad (75°). 


(18c) 


Figure 12 below shows the planar antenna created in Microwave Studio using the above 
parameters. The feeds of the two spiral arms were designed to be 2 mm apart. 



Figure 12. Illustration of a two-arm planar equiangular spiral antenna. 

With these design parameters, the widths of the arms have the same dimension as 
the spacing between the arms. This structure is self-complementary as defined by Babi- 
net’s principle [32]. Two structures are self-complementary if one is obtained from the 
other by exchanging the open and the active regions of the plane. 

The input impedance of any plane sheet antenna whose shape is the same as the 
shape of its complement (except for a trivial change of coordinates) is independent of 
frequency and is given by 





(19) 


where <^/2 is equal to 60;r in free space wl88.5Q [27]. 
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The constant impedance of a self-complementary antenna follows from the rela¬ 
tionship 

Z 1 Z 2 = (60;r)' (20) 

where Zj is the impedance of the antenna, and Z 2 is the impedance of its complementary 
structure. 

The input impedance of the spiral model (Z.„ =R + jX) is almost resistive, with 

(Z.^ with R varying between approximately 120 Q and 160 Q. Figures 13 and 

Figure 14 show the predicted variation of the input impedance of the spiral antenna over 
the frequency range from 1 GHz to 6 GHz. The difference between theory and simulation 
is attributed to the finite arm length and the finite antenna plate thickness. Notice that the 
reactance is small compared to the resistance, allowing the antenna to be well matched to 
the feeding transmission line over the designed operating bandwidth. 


Real Part of Z Matrix Coefficients in V/A 



Figure 13. Plot showing the variation of input resistance with frequency predicted 

using Microwave Studio. 
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Imaginary Part of Z Mab ix Coefficients in V/A 



Figure 14. Plot showing the variation of input reaetance with frequency predicted 

using Microwave Studio. 

Dyson [29] also discovered that consistently good radiation patterns can be ob¬ 
tained with spirals of only 1-1/4 or 1-1/2 turns, and the antenna pattern has a broad lobe 
perpendicular to the plane of the antenna. A spiral antenna was designed with two equal 
dimension arms, each with 1.5 turns. The ends of the two spiral arms were truncated to 
produce the smallest physical antenna for a given lower resonant frequency. Alterna¬ 
tively, the end of the spiral arms can be tapered, which will result in a more constant in¬ 
put impedance for the antenna [33]. Figure 15 shows the three-dimensional (3-D) radia¬ 
tion pattern of the single planar antenna created. Figures 16 and 17 show the directivity of 
the antenna for (zi = 0° and (j) = 90°, respectively. 
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Figure 15. Plot showing the 3-D radiation pattern of the spiral antenna at 1 GFlz 

predicted by Microwave Studio. 



Frequency = 1 


Tlieta / Degree 


Figure 16. Plot showing the directivity of the antenna for ^ = 0° as predicted by 

Microwave Studio. 
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Farfield ff_01.0000 [1] Directivity_Abs in dBi 



Figure 17. Plot showing the direetivity of the antenna for (j) = 90° as predicted by 

Microwave Studio. 

The radiation pattern is broad and bi-directional with equal beams radiated from 
the front and back of the structure, perpendicular to the plane of the antenna. The beam is 
circularly polarized on its axis over the usable bandwidth with the axial ratio rising 
slowly as the off-axis angle increases. Polarization is more sensitive to changes in fre¬ 
quency than is the pattern shape [37]. There is also no tilt to the lobe of the symmetrical 
antenna. By placing an absorbing material in a cavity behind the spiral antenna, the an¬ 
tenna exhibits a unidirectional radiation pattern. The disadvantage is that only half of the 
input power is transformed into radiated power because of the presence of the absorber. 
Morgan [34] considered the problem of achieving forward radiation from equiangular 
spiral antenna by rear loading. 

Since the active region rotates at a rate dependent upon the spiral parameter as the 
frequency decreases, the radiation pattern would also rotate as the frequency decreases at 
the same rate. Hence, except for rotation, the radiation characteristics are ideally inde- 
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pendent of frequeney [29]. The amount of rotation, (in radians), is given by the fol¬ 
lowing equation 


= — In 
a 


/i+A/ 

V /i y 


( 21 ) 


where a = 0.5, A/" = 0.5 GHz, giving A^ = 46°. The amount of rotation is similar to the 
simulation as shown in Figures 18 and 19. 
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Figure 18. Plot showing the 3-D radiation pattern of the spiral antenna at 2 GHz. 
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Figure 19. Plot showing the 3-D radiation pattern of the spiral antenna at 2.5 GFIz. 

The commonly used criterion for circular polarization requires that the axial ratio 
of the antenna be less than 6 dB [35]. The antenna axial ratio |AR| is defined as the ratio 
major axis electric field component to that along the minor axis. The sign of the axial ra¬ 
tio is positive for right-hand sense and negative for left-hand sense. Axial ratio is com¬ 
monly expressed in decibels as 201og|AR|. 

The idealized polarization of a right-handed spiral antenna for signal transmission 
and reception is shown in Figure 20. When used as a transmitting antenna, the radiation 
from the planar spiral antenna is right-hand circularly polarized (RHCP) in the forward 
direction (above) and left-hand circularly polarized (LFICP) in the backward direction 
(below). When used as a receiving antenna, the right-hand spiral antenna receives the 
right hand circularly polarized component (co-polarized component) of the electric field 
incident from the forward direction, while the left-hand circularly polarized component 
(the cross-polarized component) passes by the antenna mainly undetected. This circularly 
polarized field of the antenna can be beneficial when the antenna is used over the ground. 
Figure 21 shows the ratio of the electric fields of the LFICP and the RHCP components as 
a function of the angle of incidence [26]. Results are shown for dry, medium moist and 
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wet soils with the radar operating at frequency of 1.8 GHz. It is evident from the plot that 
the electric field of the reflected wave is predominantly cross-polarized to that of the in¬ 
cident wave, as long as the incident angle is not too large. Therefore, an ideal spiral an¬ 
tenna that radiates a LHCP wave towards the air/ground interface will not receive the 
wave reflected from the interface, because it is predominantly RHCP. 



Transmission 


Maximum Reception Null Reception 

for LHCP for RHCP 


Figure 20. Idealized polarization characteristics of the right-handed spiral antenna. 
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Figure 21. Ratio of £■-fields of LFICP and RFICP waves as a function of the angle of 

incidence (From Ref [26].). 

The polarization of the radiated wave is determined by the length of the arms of 
the spiral antennas. Polarization is linear when the antenna arms are very short compared 
to one wavelength, and becomes elliptical as the frequency increases. For frequencies at 
which the arm lengths of the spiral antenna are greater than one wavelength, or slightly 
less, the polarization is circular and the input impedance remains almost constant as fre¬ 
quency varies [29, 35]. A typical variation in axial ratio of the on-axis field as a function 
of frequency for a one-turn slot antenna is shown in Figure 22 [29]. 
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Figure 22. On-axis polarizations for one-turn spiral slot against frequeney 

(From Ref [29].). 


The lowest frequeney of operation is where the electrieal length of the arm is one 
wavelength. The antenna simulated here has an arm length of 0.52 m, whieh results in an 
operating frequency of 0.58 GFlz. The relationship between arm length and operating fre¬ 
quency of the spiral antenna is determined as follows. 

The total length (T) of the spiral arm can be calculated by 


i = l 



d(f) 

dp^ 


1/2 


+ 1 



( 22 ) 


^ = (A-Po) 

where and refer to the inner and outer radii of the spiral antenna, and a = 0.5 rad'^ is 
the growth rate of the spiral. Substituting/?^ =0.0021 m and/?j =0.233767 m gives the 


i+A 

a ’ 


(23) 
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spiral antenna arm length as 0.52 m, whieh results in the lowest operating frequeney of 
0.58 GHz. Figure 22 shows the variation of the antenna on-axis polarization with fre- 
queney for the spiral antenna designed. 

It ean be eoneluded from the axial ratio plots (Figures 23 to 27) that the spiral an¬ 
tenna is eireularly polarized in the 0 = 0° direetion for the frequeneies simulated. Cireu- 
lar polarization is expeeted at frequeneies up to the aeeeptable Sn parameter. At this 
limit, the size of the feed region beeomes small eompared to the feed input. 

Figure 28 shows the return loss (the magnitude of the Syz-parameter in dB) of the 
planar spiral antenna in the frequeney range from 1 to 6 GHz. The return loss was found 
to be better than -10 dB from 1 to about 5 GHz. The Smith ehart normalized to the input 
impedanee of 140 Q is also shown in Figure 29. 



FfRquRniv I 


Theta De^ec 


Figure 23. Axial ratio of the spiral antenna at frequeney /= 1 GHz. 
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Figure 24. Axial ratio of the spiral antenna at frequency/= 1.5 GFlz. 



Frequency = 2 


Tlieta / Degree 


Figure 25. Axial ratio of the spiral antenna at frequency /= 2 GFlz. 
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Farfield ff_02.5000 [1] Axial Ratio in dB 



Frequency = 2.5 


Theta / Degree 


Figure 26. Axial ratio of the spiral antenna at frequency/= 2.5 GFlz. 



Frequency = 4 


Tlieta / Degree 


Figure 27. Axial ratio of the spiral antenna at frequency/= 4 GFlz. 
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Figure 28. Return loss for planar spiral antenna. 
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Figure 29. Smith ehart plot of planar spiral antenna impedanee normalized to 140 Q. 
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Thaysen [33] has also observed that there is very little difference in the Su- 
parameter for the spiral antenna with or without a cavity. The bandwidth of the antennas 
is often defined by the range of the frequencies that mark the bounds of acceptable per¬ 
formance. The above antenna design has a bandwidth of approximately 5 GHz. 

The equiangular spiral antenna described is a balanced antenna and requires a 
balun to transform the unbalanced coplanar waveguide (CPW) feed line to a balanced 
coplanar strip (CPS) feed line for input to the spiral antenna. The bandwidth of the sys¬ 
tem is limited by the balun when an electrically unbalanced feed line is used [37]. Thay¬ 
sen [38] presented a design for a wideband balun to accomplish the transition from the 
unbalanced CPW to the balanced CPS transmission line. The balun has insertion loss of 
less than 3 dB from 100 kHz and up to a frequency of 3.85 GHz. A four-section Cheby- 
shev impedance transformer was designed in order to transform the impedance from 50 
ohms to the spiral antenna impedance of 80 ohms. The reflection coefficient was found to 
be less than 5 dB. 

Simulations of spiral antennas placed on substrates of different relative permittiv¬ 
ity have shown that the impedance of the antenna decreases with increased permittivity. 
Figure 30 shows that the input impedance of the antenna for a substrate permittivity of 
1.06 is 200 Q at 2 GHz. 



Frcqucncr.' Qk 


Figure 30. Variation of input resistance with frequency for =1.06. 
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c. 


SUMMARY 


A frequency-independent antenna is one whose radiation pattern, impedance and 
polarization remain virtually unchanged over a large bandwidth. The design of a fre¬ 
quency-independent antenna is based upon the fundamental principle that, if the shape of 
the antenna were such that it could be specified entirely by angles, its performance would 
be independent of frequency. 

The planar spiral antennas have useful properties for subsurface applications. The 
polarization of the radiated wave from the planar spiral is determined by the length of the 
arms of the spiral antennas. For frequencies where the arm lengths of the spiral antenna 
are greater than one wavelength, the polarization is circular and the input impedance re¬ 
mains almost constant as frequency varies. The impedance of the antenna decreases with 
increased substrate permittivity. 

The counter-wound spiral antenna (CWSA) design examined in the next chapter 
aims to minimize losses due to polarization mismatch between the scattered fields from 
the subsurface and the receiving antenna. When the scattered field from the buried ob¬ 
jects or interface is linearly polarized, the use of a circularly polarized receiving antenna 
would result in 3-dB polarization loss. 
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IV. COUNTER-WOUND SPIRAL ANTENNA DESIGN 


As mentioned in Chapter III, the ideal equiangular spiral antenna will not receive 
waves reflected from the air/ground interface because the electric field (£'-lield) of the 
reflected wave is predominantly cross-polarized to that of the incident wave, as long as 
the incident angle is not too large. 

The RCS of a target depends on many factors such as the frequency and polariza¬ 
tion of the incident wave, as well as the target material and target aspect [39]. Target as¬ 
pect refers to the orientation of the target with respect to the £’-field of the incident wave. 

The polarization of a uniform plane wave describes the shape and locus of the tip 
of the £'-field vector in the plane orthogonal to the direction of propagation, at a given 
point in space as a function of time. Figure 31 [40] shows a snapshot of a plane wave 
propagating in the +y-direction at a fixed instant of time. The wave is linearly polarized 
since the ^-field at a fixed point oscillates up and down along a vertical line in the z- 
direction. 



Figure 31. Spatial variation of E-field at t = 0 for the plane wave propagating in 

+y-direction (From Ref. [40].). 
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For a plane wave propagating in the z-direetion, the z-eomponents of both the 
eleetric and magnetie fields are zero. The E-field of a time-harmonic wave may be de¬ 
scribed by phasors (denoted by ~) 

E{z) = xE^ (z) + yE^ (z) (24) 

with 

E,(z) = (25a) 

E,(z) = E,„e-i‘-e‘\ (25b) 

where E^q and E^^mq the amplitudes of the components of E^{z) and E^{z), respec¬ 
tively. The phase of a wave is defined relative to a reference condition with specified val¬ 
ues of z and t. The wave polarization depends on the relative phase of E^with respect 

to Ey, not their absolute phases. 

For a linearly polarized wave, the E-field at a fixed point oscillates back and forth 
along a line, that is, the E^{z,t) andE^(z,t) components are either in phase or out of 
phase {S = 0,7r). For a circularly polarized wave, the components of the E-fields, 
E^ (z, t) and E^ (z, t) remain constant and the phase difference between the fields is 90 

degrees {5 = ±;r/2). Hence, a circularly polarized wave can be obtained by superimpos¬ 
ing two equal-amplitude linearly polarized plane waves that are in space and time quadra¬ 
ture. 

The sense of the polarization is defined in terms of the rotation of the E-field as a 
function of the time in a fixed plane orthogonal to the direction of propagation, which is 
opposite of the direction of rotation of E-field as a function of distance at a fixed point in 
time. The designation of RHCP is determined by the right-hand rule and FHCP by the 
left-hand rule. Figure 32 shows RHCP and FHCP waves. 
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Figure 32. F’-field variations of the RFICP and the LFICP waves (From Ref [40].). 

Polarization effieieney is given by the polarization mismatch factor, p. It is used 
to determine the polarization mismatch between the received wave and the antenna. The 
polarization efficiency varies from a minimum of 0 to a maximum of 1 as the incoming 
wave and the receiving antenna vary from a completely mismatched to a completely 
matched condition. A complete match {p = \) exists when the wave and the antenna po¬ 
larization states are exactly identical, which will result in maximum antenna output volt¬ 
age. A complete mismatch of the system (;? = 0 ) occurs when the wave and antenna are 
cross-polarized, which will result in a null output. Orthogonal linear states such as hori¬ 
zontal and vertical linear polarization (LP) and right-hand and left-hand circular polariza¬ 
tions (CP) are examples of cross-polarized states. 

The polarization of the scattered field is not necessarily the same as that of the in¬ 
cident field. Most targets will generate cross-polarized scattering components due to mul¬ 
tiple reflections and diffractions arising from the nature of the objects. A scattering ma¬ 
trix can be used to completely specify the relationship between the incident and the scat¬ 
tered fields 
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where {Sp^} are the seattering parameters, and p and q ean be either 6 or (j). The first 

parameter speeifies the polarization of the reeeive antenna and the seeond the polarization 
of the wave. The elements of the seattering matrix are eomplex quantities. 

One half of the power available from a LP wave is lost when reeeived by a CP an¬ 
tenna. The LP wave ean be deeomposed into equal eomponents of RHCP and LHCP 
waves. If the reeeiving antenna is eireularly polarized, it will only be reeeiving either one 
of the components, since the other one will be co-polarized. The same is true for a CP 
wave and an LP antenna. In most applications, this 3-dB loss is significant and hence a 
receiving antenna that is matched to the polarization of the incoming wave must be used. 

Thus, for an antenna to be able to detect a buried object, the scattered electric 
field from the object must contain a significant co-polarized component. For a buried per¬ 
fectly conducting thin rod shown in Figure 33, the scattered field from a thin rod will be 
linearly polarized, and its field can be decomposed into equal LHCP and RHCP compo¬ 
nents. For the spiral antenna shown, only half of the backscattered power will be received 
by the antenna, reducing the probability of target detection. 




LHCP 


LHCP RHCP 



Figure 33. Geometry of a monostatic GPR that uses a single spiral antenna to detect 

thin rods buried in the ground. 
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A, 


COUNTER-WOUND SPIRAL ANTENNAS IN FREE SPACE 


The counter-wound spiral antenna (CWSA) design examined here aims to mini¬ 
mize losses due to polarization mismatch between the scattered fields from the subsurface 
and the receiving antenna. This is achieved by having two sets of dual-arm spiral anten¬ 
nas, counter-wound and stacked over each other along the antenna main axis as shown in 
Figure 34. The design of each of the planar duo-arm equiangular-spiral anteimas was 
based on the same fundamental design principles discussed in Chapter III. The two planar 
anteimas can be fabricated on separate sides of the same substrate, with the substrate 
thickness determining the spacing between the two spirals. 



Figure 34. Geometry of a counter-wound spiral antenna. 

The objective of the design is to achieve an antenna radiation pattern and imped¬ 
ance that remains virtually unchanged over a large bandwidth. With the CWSA, the po¬ 
larization of the wave transmitted by the antenna can be made linear, either vertically or 
horizontally polarized, by varying the delays to the signals that are being sent to the an¬ 
tenna feeds of the spirals. This design allows the buried objects to be radiated by both 
vertically and horizontally polarized waves, reducing the polarization loss in the received 
scattered signal from the subsurface objects, thereby allowing better detection of targets 
which may lie in any arbitrary orientation with respect to the main antenna axis. 
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Although the linearly polarized wave will reeeive refleetions from the ground, 
they ean be removed by signal proeessing, using delay line eaneellers or other teehniques 
eommonly used in today’s airborne radar systems. Peter [41] has sueeessfully aehieved 
greater than 30 dB of elutter suppression by subtraeting the extraneous elutter (without 
the target) eoherently from the data colleeted with the desired target present. This signal- 
to-elutter ratio (SCR) was further enhaneed by using high range resolution proeessing. 

The performanee of several CWSAs was determined through simulation for dif¬ 
ferent separations and orientations between the two planar spiral antennas. The resultant 
polarization of the CWSA was also determined by ealeulating the axial ratio of the E- 
fields at different distanees and displaeements from the antenna axis. The frequency- 
independent characteristics of the antenna, like its impedance and radiation patterns over 
a wide bandwidth, are also determined through software simulation. The results of the 
simulations for different parameters in the antenna design are shown in the following 
subsections. 

When both antennas are transmitting, the combined signal will be linearly polar¬ 
ized under ideal conditions, that is, when the two spirals lie very close together and yet do 
not interfere with each other. Similarly, on receiving, the combined output of the two spi¬ 
rals will be linearly polarized. 

As shown in Figure 35, the top spiral arms, which have a bi-directional radiation 
pattern, will transmit a LHCP wave to the rear counter-wound spiral antenna. Ideally the 
wave should pass through the rear antenna because it has the opposite polarization sense. 
The addition of the E’-fields from the top antenna (RHCP) and the rear antenna (LHCP) 
will result in a linearly polarized wave in the forward direction. 
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Figure 35. Setup for determining the reflection coefficients of the counter-wound 

spiral antenna. 

Figure 36 shows the 5'/;-parameter of the top planar spiral antenna for a substrate 
thickness of 20 mm, for frequencies ranging from 1 to 6 GHz, with the rear antenna feed 
terminated in a resistive load equal to the estimated antenna impedance of 140 Q. This 
condition gives an indication of the coupling or interference presented by the bottom an¬ 
tenna to the top antenna. The 6';/-parameter of the antenna was found to be better than 
-10 dB from 1 GHz to 3.8 GHz. The subsurface radar can then be designed to operate in 
this frequency band for low reflection coefficients. The Smith chart normalized to an im¬ 
pedance of 140 Q is also shown in Figure 37. For this setup, the top and rear spiral arms 
are rotated at 90 degrees as shown in Figure 34. 

It can be seen from these plots that an operating bandwidth of about 3.8 GHz is 
possible for an antenna separation of 20 mm. Figures 38 and 39 shows the simulation re¬ 
sults for a substrate thickness of 30 mm. It is concluded from the plots that a substrate 
thickness of 30 mm is required to achieve a continuous bandwidth of about 4.8 GHz 
(based on Su alone) when the counter-wound spiral antennas are setup as in Figure 34. 
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Figure 36. Return loss of the CWSA with substrate thickness of 20 mm. 



Figure 37. Smith chart for the CWSA with substrate thickness of 20 mm. 
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Figure 38. Return loss of the CWSA with substrate thickness of 30 mm. 



Figure 39. Smith chart for the CWSA with substrate thickness of 30 mm. 
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Simulation has also been performed to determine the performanee of the CWSA 
design as shown in Figure 40. In this design, the spiral arms are eounter-wound, but are 
not rotated 90 degrees to each other as in the previous case. Simulations were performed 
to determine the return loss of the un-rotated CWSA for a substrate thickness of 100 mm, 
for frequencies ranging from 1 to 6 GHz. The rear antenna feed was also terminated with 
a resistive load value equal to the estimated impedance of the spiral antenna. 

The return loss and Smith chart of the un-rotated CWSA are as shown in Figures 
41 and 42. It can be concluded from the plots that a substantial increase in the separations 
between the planar spirals is required in order to achieve similar performance to that of 
the rotated case. This un-rotated design allows a smaller overall antenna footprint but at 
the expense of a much wider separation between the spiral arms. 



1 ii/i'r lypr a PFC 


Figure 40. Geometry of an un-rotated counter-wound spiral antenna. 
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Figure 41. Return loss of the un-rotated CWSA with substrate thickness of 100 mm. 
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Figure 42. Smith chart for the un-rotated CWSA with substrate thickness of 100 mm. 
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The polarization when both planar duo spiral arms are simultaneously excited is 
shown in Figure 43. The £’-field of the RHCP wave radiated by the front planar spiral is 
being added to the £’-lield of the LHCP wave transmitted by the counter-wound rear spi¬ 
ral. The resultant polarization of the CWSA would be a linearly polarized wave traveling 
in the forward direction towards the ground. 



Figure 43. Linear polarization from a counter-wound spiral antenna. 


The linear polarization can be shown mathematically as follows. The £’-lield of a 
RHCP wave can be represented as a vector given by 


ERHCP=Eoix-jy)e 


-Az-f 


(27) 


and the £’-lield of a LHCP wave represented by 


-jk[z+E\ 

ELHCP=Eoix + jy)e ^ "V 


(28) 


where x and y are the unit vectors in the x and y directions respectively. The separation 

between the front and rear spiral antennas is given by d, and is an arbitrary phase dif¬ 
ference between the front and rear antennas. This is illustrated in Figure 43. 
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For an observation point at z, the resultant £’-field due to the addition of the RHCP 
and LHCP waves is given by 
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(29) 


whieh is a linearly polarized wave with its eleetrie field veetor rotated from the x-axis by 
an angle a where 


tan a = tan 


kd (b 
— + — 
2 Ij 


(30) 


A phase differenee (j) = 180° between the front and rear antennas will result in a rotation 
of the resultant field Ej,qj by 90°. For the ease of J = 0 and = 0°, 


Ejqj = 2E^^e , 

whieh is a horizontally polarized wave. For the ease of J = 0 and (j) = 180° , 


(31) 


^TOT 


= 2E^e 


-jkz 




(32) 


whieh is a vertieally polarized wave. 

In order to generate a horizontally polarized wave, the phase of the front spiral’s 
exeitation signal was delayed by an amount proportional to the wavelength of operation 
and the separation between the two planar antennas. The relationship between the phase 
delays (^ ), wavelength (A) and separation (d) is given by 


(j) = 


2nd 


rad. 


(33) 
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For a frequency of operation of 1.2 GHz and separation of 10 mm, the required 
phase delays to be applied to the front planar antenna in order to achieve a horizontally 
polarized wave would be -14.4°. Figure 44 shows the amount of phase delay required 
for the same antenna separation operating at different frequencies. 
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Figure 44. Phase delays required for an antenna separation of 10 mm. 

The polarization of the fields at the antenna axis was verified through simulation. 
Figure 45 shows the vector plot of the resultant ^-field components on the antenna plane 
at a distance of about 19.3 mm away from the antenna. Figures 46 to 48 show the indi¬ 
vidual ^-field components in the x, y and z-directions. The antenna was found to be line¬ 
arly polarized near the antenna axis. For the same delay settings, the polarization orienta¬ 
tion was verified at a frequency of 2.4 GHz. Figures 49 to 52 show the ^-field plots at 
that frequency with the front planar antenna delayed hy (j) = -14.4°. The polarization of 
the antenna was also found to be linear at a frequency of 2.4 GHz. The far-field axial ra¬ 
tio plots for frequencies from 1.9 to 2.4 GHz are shown in Figures 53 to 58. 
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Figure 45. Vector plot of the resultant E'-fields at frequency /=! .2 GHz. 
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Figure 46. ^-field component in the x-direction at frequency/=1.2 GHz. 
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Figure 47. ^-field component in the j-direction at frequency/ =1.2 GHz. 
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Figure 48. ^-field component in the z-direction at frequency/=! .2 GHz. 
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Figure 49. 


Figure 50. 
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Figure 51. ^-field component in the j-direction at frequency/= 2.4 GHz. 
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Figure 52. ^-field component in the z-direction at frequency/= 2.4 GHz. 
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Figure 53. Axial ratio of the antenna at frequency/= 1.9 GHz. 
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Figure 54. Axial ratio of the antenna at frequency/= 2 GHz. 
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Faifield ff_02.1000 [1(1.0,-14]+2[1.0^.0],[3]] Axial Ratio in dB 
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Figure 55. Axial ratio of the antenna at frequeney / = 2.1 GHz. 
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Figure 56. Axial ratio of the antenna at frequeney / = 2.2 GHz. 
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Faifield ff_02.3000 [1(1.0,-14]+2[1.0^.0],[3]] Axial Ratio in dB 
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Figure 57. Axial ratio of the antenna at frequency/= 2.3 GHz. 
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Figure 58. Axial ratio of the antenna at frequency/ = 2.4 GHz. 
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As described previously, the polarization can be changed from a vertically polar¬ 
ized wave to a horizontally polarized wave by adding an additional phase delay of 
180° degrees to the front spiral of the CWSA. A simulation was performed by changing 
the phase delays of the front planar antenna to -104° . The ^-field components are shown 
in Figures 59 to 62. The 90° phase shift rotated the plane of polarization by 90° as ex¬ 
pected. 
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Figure 59. Vector plot of E-fields at frequency of 2.4 GHz with delay of -104° . 
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Figure 60. ^-field component in the x-direction at frequency of 2.4 GHz with -104° 

delay. 
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Figure 61. ^-field component in the y-direction at frequency of 2.4 GHz with -104° 

delay. 
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Figure 62. ^-field component in the z-direction at frequency of 2.4 GHz with -104° 

delay. 


B, SUMMARY 

Several important aspects of the performance of the CWSA have been verified 
through a series of simulations. It has been shown that the rotated CWSA (with the plane 
of the front and rear spiral antenna planes rotated 90° to each other along the main an¬ 
tenna axis) gives a better performance in terms of reflections as compared to the un¬ 
rotated case. The CWSA far-field polarization has also been verified to be linearly polar¬ 
ized on axis for any amount of delay applied to the input of the front planar antenna. The 
tilt angle of the linearly polarized wave can be changed from vertically polarized to a 
horizontally polarized by increasing the signal excitation delay to the front antenna by 
180°. This allows buried objects with unknown aspect angle with respect to the receiving 
antenna to be detected and identified by the radar system without large polarization loss. 

Use of the CWSA can be extended to other applications where a large antenna 
bandwidth is required. The CWSA retains the wideband antenna characteristics of the 
spiral antennas with an improved detection capability by using linearly polarized (both 
vertical and horizontal) waves. 
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V. CONCLUSIONS AND FUTURE WORK 


A. CONCLUSIONS 

This thesis has addressed many issues pertaining to the design of subsurfaee ra¬ 
dar systems and wideband antenna design. In partieular, the researeh investigated the fae- 
tors eontributing to the performanee of a subsurfaee radar system, its arehiteeture, and 
design of frequeney-independent antennas. The result is a new antenna design, the 
eounter-wound planar spiral, which is capable of variable linear polarization over a wide 
bandwidth. 

Frequency-independent antennas refer to antennas that have no theoretical limita¬ 
tion on the bandwidth of operation. In practice, they are antennas whose radiation pattern, 
impedance, and polarization that remain virtually unchanged over a large bandwidth. A 
spiral antenna offers the advantages of having a large operating bandwidth, constant real 
impedance and a circularly polarized wave, which automatically reduces the amount of 
reflection received from the air/ground interface due to cross-polarization. The signal re¬ 
ception is, however, constrained by the target orientation with respect to the antenna axis. 
Only the co-polarized portion of the reflected wave from the object will be received by 
the antenna, which may result in 3-dB (6-dB for both transmit and receive) polarization 
loss, which is not acceptable in many applications. 

This thesis examined the design and performance of a counter-wound spiral an¬ 
tenna. The CWSA can produce a linearly polarized wave in the forward direction by 
varying the phase of the excitation to the front planar spiral arms. The phase delay is de¬ 
pendent on the separation between the planar spiral arms and the frequency of operation. 
The resultant wave from the antenna will be linearly polarized. The plane of polarization 
can be made either horizontal or vertical by varying the delay to the antenna feeds. This 
design allows buried objects with unknown aspect angle with respect to the receiving an¬ 
tenna to be detected by the radar system without large polarization loss. Variable linear 
polarization allows the extraction of the maximum amount of information from the scat¬ 
tered echo. The reflections due to the returns from the air/ground interface can be reduced 
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by using signal processing techniques. The CWSA design may be extended to other ap¬ 
plications where large antenna bandwidth is required. 


B, FUTURE WORK 

The counter-wound spiral antenna can be fabricated on a substrate, with its per¬ 
formance measured in the anechoic chamber. The antenna patterns of the fabricated an¬ 
tenna can then be measured and compared with the theoretical values and simulation re¬ 
sults. The impedance match between the antenna feeds and the spiral antenna could be 
measured with test equipment and its effect on the overall antenna performance deter¬ 
mined through measurements. However, it will be necessary to design a broadband balun. 

For the model discussed here, the air/ground interface was a smooth, planar sur¬ 
face. In practical applications, the interface could be rough and may even support various 
forms of vegetation. The effects of these factors on the performance of the monostatic 
GPR using CWSA can then be determined by measurement. The benefits of using this 
antenna structure for bistatic GPR have not been examined and is worthy of investigation. 
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